Feshbach Resonances in Cesium at Ultra-low Static Magnetic Fields 
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We have observed Feshbach resonances for 133 Cs atoms in two different hyperfine states at ultra- 
low static magnetic fields by using an atomic fountain clock. The extreme sensitivity of our setup 
allows for high signal-to-noise-ratio observations at densities of only 2 x 10 7 cm -3 . We have repro- 
duced these resonances using coupled-channels calculations which are in excellent agreement with 
our measurements. We justify that these are s-wave resonances involving weakly-bound states of 
the triplet molecular Hamiltonian, identify the resonant closed channels, and explain the observed 
multi-peak structure. We also describe a model which precisely accounts for the collisional processes 
in the fountain and which explains the asymmetric shape of the observed Feshbach resonances in 
the regime where the kinetic energy dominates over the coupling strength. 

PACS numbers: 67.85.-d,34.50.Cx,37.10.Vz,06.30.Ft 



The achievement of Bose-Einstein condensation |l| has 
stimulated remarkable developments in atomic physics. 
Ultracold atoms have found applications in metrology [2| 
and high-precision measurements of physical constants 
Q; they can be cooled down to quantum degeneracy 
and used to simulate condensed-matter systems 0, Hj . 
A fundamental feature of ultracold atomic gases, un- 
derlying most of their present applications, is that the 
interparticle interactions can be tailored at will, using 
scattering resonances that occur in low-energy collisions 
between two atoms 0|. These Feshbach resonances are 
usually obtained using an external static magnetic field 
0. Their accurate characterization is intimately linked 
to a detailed knowledge of the interatomic interaction [H 
and involves coupled-channels calculations [9j. 

We report on the measurement of multiple Feshbach 
resonances in 133 Cs using an atomic fountain clock, 
and present their theoretical characterization using the 
coupled-channels method. The extreme accuracy of fre- 
quency measurements in modern atomic clocks provides 
the means to reveal effects of atomic collisions in a regime 
of very weak interactions. The excellent agreement be- 
tween experimental measurements and theory confirms 
that the interaction between Cs atoms is now well un- 
derstood and modeled. The resonances that we analyze 
are unusual for two main reasons. First, they occur at 
magnetic fields of the order of a few milliGauss, which 
makes them the lowest-static-field resonances investi- 
gated up to now. In these ultralow magnetic fields, the 
quasi-degeneracy of all collisional channels with a triplet 
two-atom electronic spin plays a key role and conveys a 
multi-peak structure to the resonances. Second, we have 
measured them in a regime where the kinetic energy dom- 
inates over the resonance width. In this regime, they ap- 
pear in the magnetic field dependence of the clock shift 
as asymmetric features which occur close to the zero- 



temperature resonant field. 

The further experimental characterization of these low- 
field resonances using density-independent interferome- 
try [lQj |. combined with the enhanced sensitivity to the 
values of fundamental constants near a Feshbach reso- 
ld, could be used to probe the constancy of 
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the proton-to-electron mass ratio and the fine structure 
constant. Furthermore, these resonances involve atoms in 
two different spin states and thus pave the way towards 
the study of quantum magnetism in ultracold Cesium 
gases containing two different hyperfine states. 



-4 m f =0 +4 
f= 4 _._ 

f=3 J 1 

-3 m f =0 +3 

(2) 

J (3) 



FIG. 1. Left: schematics of the atomic fountain. Right: 
hyperfine levels used in the experiment. Populations after 
launch and transitions excited for the state selection (1). A 
horizontal laser beam then pushes away / = 4 atoms. Popula- 
tions at the start (2) and during (3) the Ramsey interrogation. 

Experimental setup. — The experiment is done in a 
fountain geometry which has already between described 
extensively (see e.g. [2[) and which is sketched in Fig. [U 
The frequency of the \ f — 3,ra/ = 0) — > \ f = 4, ra/ 
0) hyperfine transition is probed during the ballistic flight 
of a cloud of 133 Cs atoms laser-cooled to ~ 1 /iK. The 
Ramsey interrogation occurs between the upward and 
downward traversals of a microwave cavity. After the 
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Ramsey interrogation, the atom numbers in the / = 3 
and / = 4 hyperfine states are measured by laser-induced 
fluorescence detection in order to determine the transi- 
tion probability. Before the interrogation, state selection 
is applied to the up-going cloud by means of microwave 
and laser interactions (see Fig. [1]). For the present ex- 
periments, we select not only the |3,m/ = 0) clock state 
but also an additional |3, rrif ^ 0) state, and measure the 
clock frequency shifts due to this state. 

Clock shift measurements of Feshbach resonances — 
Collision-induced frequency shifts depend on elementary 
collisional properties but also on the atomic spatial and 
velocity distributions. This latter dependence is even 
stronger in the present experiment, first because of the 
evolution of the atomic cloud during the Ramsey inter- 
rogation (see e.g. [13]), and second because we are in a 
regime of strong sensitivity of the measured shifts to the 
collision energy. We determine collision shift ratios in 
a way that minimizes the impact of atomic distributions 
which are difficult to control with high precision. We per- 
form interleaved frequency measurements with 3 config- 
urations, leading to 3 measured frequencies: 
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and z^o-m/- Firstly, the |3,ra/ = 0) state is selected with 
the maximum possible atom number No. Secondly, the 
1 3, ra/ = 0) state is selected with the atom number No/2. 
Thirdly, No atoms in |3,m/ = 0) are selected together 
with N mf atoms in another chosen |3,m/) state, as il- 
lustrated in Fig. [I] for m/ = —3. The expanding atomic 
cloud is truncated during the Ramsey cavity tr aver sals, 
so that the detected atoms are only a fraction (~ 20%) 
of the initially selected atoms. We choose to character- 
ize the number of atoms using the detected atoms, hence 
No, No/ 2 and N mf will refer to the atom numbers as 
measured in the detection. For No and No/2, this is the 
sum of |3, 0) and |4, 0) atoms since some atoms which are 
initially in the state |3,0) are excited to the state |4,0) 
during the Ramsey interrogation. A crucial feature of 
our experiment is to perform the microwave excitation 
for state selection with the (interrupted-) adiabatic pas- 
sage method described in [14] in order to ensure quasi- 
identical space and velocity distributions for all states 
(m/ = and m/ ^ 0) and all configurations, notably 
the first and the second one. We can prepare the third 
configuration with any of the six m/ = ±1, ±2, ±3 states. 

In a given configuration, the frequency shift of the 
|3,0) — > |4,0) transition is given by: 

5v = n poKo(B,V ryV ) +n mf p mf K rnf (B,V ryV ) , (1) 

where no and 

rtffi j are the detected atom numbers, po and 
p mf are the effective densities per detected atom, and Ko 
and K mf are the collision shifts scaled to the effective 
densities. These functions include collisional properties, 
which depend on the magnetic field B. They also depend 
on the space and velocity distributions V r ^ Vl and more 
generally on the fountain geometry. Starting from the 
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measured frequency shifts and the detected atom num- 
bers, we compute the shift per detected m/ = atom, 
Ao ; o, and the additional shift A mf; o due to the m/ ^ 
population, per detected m/ ^ atom: 
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Our state selection ensures quasi-identical distributions 
for all states, so that po — p m/ , and, hence, A m/; o/Ao ; o — 
K mf (B,V riV )/Ko(B : V riV ). This quantity is as close to 
intrinsic collisional properties as possible in our experi- 
ment. Notably, it does not depend on the detected atom 
numbers No and N mf . Typically, No ~ 5 • 10 6 and 
N mf « Ao- The corresponding effective density during 
the Ramsey interrogation, Aopo ~ 2 • 10 7 cm -3 , is many 
orders of magnitudes lower than in typical quantum gas 
experiments. The mean free path is ~ 35 m and the 
mean time between collisions is ~ 5000 s, i.e. 3 orders of 
magnitude longer than the experimental cycle. 

We have determined Ao ; o, A mf; o and R mf] o = 
A mf .o/Ao-o for all m/ states as a function of the mag- 
netic field B from to 100 mG. The magnetic field is 
known via the spectroscopy of the first-order-sensitive 
|3,m/ = 1) — > |4, 777/ = 1) transition. It is stable to 
~ 40 nG and homogeneous to better than 10 -2 . Our 
measurements of R mf] o(B) are shown in Fig. [2l The 
magnetic field keeps the same downward orientation over 
the entire height of the fountain to avoid Major ana tran- 
sitions and to ensure a good control of the quantization 
axis. Under these conditions, selecting a —777/ state for 
a measurement is equivalent to probing the +777/ state 
with the field — B. Therefore, we plot measurement re- 
sults with B < which are, in fact, taken with a negative 
m/ state. For all 3 states, we observe a dramatic depen- 
dence of R mf ;0 on B. Instead, we measure no significant 
change of the clock collision shift Ao ; o(B), at a level lim- 
ited by the dependence of this quantity on the effective 
density po> which itself is quite sensitive to variations 
in the atomic distributions T> TjV . Within these limits, 
Ko(B,V r ^ v ) remains constant over the entire range of 
our experiments. It is equal to the large negative clock 
shift which affects Cs fountain clocks [15l-ll7j. Hence, the 
observed behavior of R mf; o(B) relates to Km f (B,7) riV ), 
which we attribute to Feshbach resonances either in the 
|3,0;3,777/) or the |4,0; 3,777/) channel. 

The precise control of the magnetic field and the high 
signal-to-noise ratio of the data allow for a stringent 
comparison to two theoretical approaches: (i) a coupled- 
channels calculation of the scattering length characteriz- 
ing interactions at zero temperature as a function of B, 
and (ii) a finite-temperature model of the clock collision 
shift in the fountain geometry, which explains the asym- 
metric shape of the observed resonances. 
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FIG. 2. Feshbach resonances at low static magnetic fields. Left: Mf — 3; middle: Mf — 2; right: Mf — 1. Lower graphs: 
measured collision frequency shift ratio Rm f ;0 derived from the experiment as a function of the static magnetic field. Upper 
graphs: theoretical values for the scattering length. The vertical dashed lines show the theoretical resonant field values. 
Experiment (effective temperature ~ 900 nK) and theory (T = 0) are completely independent (see text). 




Calculation of the scattering length. — We describe the 
system in the center-of-mass frame of the atom pair. 
Neglecting the spin-spin interaction, which yields no sig- 
nificant contribution to our observables, the interaction 
is spatially isotropic. We limit our analysis to s-wave 
interactions governed by the following Hamiltonian [9[: 

H=^ + Vei(r) + VM + V z , (3) 

where r is the interatomic distance, p is its conjugate 
momentum, and p = ra/2 is the reduced mass of the 
atom pair. The central part of the interaction is given 
by V e \(r) = Vs(r)Ps + Vt(^)Pt, where Ps and Pt are 
the projectors onto the electronic-singlet and triplet sub- 
spaces. The term Vhf = &hf(si • i\ + s<i • 12)/^ is the hy- 
perfine interaction, where Sj and ij are the spin operators 
of the electron and the nucleus of atom j. The operator 
Vz = S z is the Zeeman term [l8|, with fig being 

the Bohr magneton and S z = s\ z + S2 Z being the total 
electronic spin projection along the quantization axis e z . 

We calculate the magnetic field dependence of the scat- 
tering length associated with the zero-energy scatter- 
ing state corresponding to the levels populated in the 
experiment. The Hamiltonian H conserves the pro- 
jection Mp of the total two-atom spin F = f\ + 
/2, where fj = Sj + ij is the total spin of atom 
j. Therefore, this scattering state has a definite 
value of the total spin projection Mp, on which the 



scattering length aM F (B) depends. For large inter- 
atomic separations, the atoms are in the Zeeman-dressed 
state related to the (Bose-symmetrized) two-atom state 
|/i=4, mi = 0; f'2 = 3, 77i2 = Mf ), where the quantum 
numbers fj and rrij define the magnitude and projection 
of the total spin fj. The experimental results shown in 
Fig. [2] correspond to Mp = 3, 2, and 1, respectively. 

The scattering state | ^m f ,b ) has 10 coupled compo- 
nents if Mf = 3, and 13 and 14 components for Mp = 2 
and Mp = 1, respectively. We evaluate it numerically 
using the coupled-channels approach our implemen- 
tation of which is described in [l9|. The accumulated- 
phase boundary condition [9] is applied at ro = 20 ao, 
and the asymptotic behaviour of the zero-energy scat- 
tering state is enforced at r max = 1000 a . All result- 
ing differential systems are solved using Stoermer's rule 
with adaptive stepsize control [20]. The values used for 
the accumulated-phase parameters, the hyperfine inter- 
action constant ahf , and the electronic potentials Vs and 
Vt are the same as those used in (2lj . 

Our results for the 8-wave scattering length aM F {B) 
are shown in Fig. [3j for Mp = 3, 2, and 1. The occur- 
rence of inelastic processes (such as the decay towards the 
lower-energy states having f 1 =f 2 = 3) causes a to have 
a non-vanishing imaginary part [22| and the resonances 
appear as smooth dispersive features (rather than as the 
divergences of the lossless case). The calculated reso- 
nance positions, corresponding to the minima of Im(a), 
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TABLE I. Measured and calculated resonance positions. 



are shown in Table HI The calculated positions for the 
broadest resonances compare favorably to those deter- 
mined from the experimental clock-shift measurements 
(Fig. [2]). The predicted multiple-peak structure is clearly 
visible in the experimental data for Mp = 2. 

Our numerical analysis includes only s-wave interac- 
tions, and the fact that it recovers the measured reso- 
nance positions proves that these are s-wave resonances. 
The triplet potential Vp supports a very weakly bound 
state, with the binding energy \Ep\ = ft 2 /(2fia T ) « 
ft -5 kHz = /iB-4mG, where ap = 2400 ao is the scattering 
length associated with Vp |8|. For a given value of Mp, 
the two-atom internal states | fa = 4, fa = 3, F, Mp ) are 
electronic-triplet for all allowed odd values of F. For B = 
0, each of these triplet channels supports the weakly- 
bound triplet state, yielding degenerate bound states 
(energy — |£r|), where is the number of triplet chan- 
nels with the quantum number Mp. For non-zero, albeit 
small, magnetic fields, the coupling due to Vz lifts this 
degeneracy, and these states cross the threshold for 
different values of B, causing multiple resonances. For 
Mp = 3 or 2, there are = 3 triplet channels (F = 7, 
5, or 3), which correspond to the three predicted reso- 
nances in these two cases. For Mp = 1, there are N± = 4 
triplet states (F = 7, 5, 3, 1); however, our coupled- 
channels results only show three resonances, probably 
because the fourth one is too narrow to be resolved in 
the presence of the three other peaks. This multiple- 
resonance physics only occurs for small B: indeed, for 
values of B larger than a few \Et\//j>Bi the Zeeman term 
Vz causes the bare weakly-bound triplet states to dis- 
solve into the continuum. 

Feshbach resonances in a fountain geometry. — To clar- 
ify the impact of finite temperatures, the atomic spatial 
and velocity distributions V r , v , and the fountain geom- 
etry, we have evaluated the clock shift using a simple 
model for the ^-matrix elements S aj (k) and S^(k) de- 
scribing the interaction between the clock states, a = 
|3, 0) and j3 = |4, 0), and the additional state 7 = |3, rrif). 
The elementary clock shift due to 7 is : 
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with p 1 being the local density, both in time and space, of 
atoms in the state 7, and k = p/% being the wavevector 



for the relative motion of the two colliding atoms. For a 
resonance occurring in the ^7 channel, we take S ai {k) = 
1 and assume that Sp 7 (k) is given by: 



E - v{B) + ±ir e 



(5) 



where E = ft 2 /c 2 /2/i is the relative kinetic energy of the 
colliding pair, v(B) is the energy detuning from the reso- 
nance, and T e = kC e is the elastic width of the resonance 
[23|| , the coupling strength C e being constant. We have 
omitted the inelastic contribution to the width, zI\/2, 
in the denominator of Eq. (|5j), as our coupled-channels 
results imply that Ti/C e <C k. The total clock shift 
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FIG. 4. Left: clock shift as a function of B, numerically 
calculated using Eqs. |4] and El for various coupling strengths 
C e - The value C e = Co = AE vec /k vec (red) is close to the 
experimental situation; results for C e = 100 Co (black), C e = 
10 Co (magenta), and C e = 0.1 Co (green) are also shown. 
Inset: close-up of the red curve for small values of B. Right: 
fit of our model (Eqs. [4] and EJ to the Mf = 3 measurements. 

is obtained by averaging Eq. (jU) over the atomic space 
and velocity distribution V rjV measured in the experi- 
ment. We calculate it using a Monte-Carlo simulation 
accounting for the collisional energy distribution (cor- 
responding to the effective temperature ~ 900 nK), the 
decrease of the atomic density with time, and the trun- 
cation of the atomic cloud in the microwave resonator. 
Figure [4] (left) shows the total clock shift as a function 
of B for various coupling strengths C e . As an exam- 
ple, we take v = 2\ib{B — Bo), corresponding to a scat- 
tering length a(B) = -C e /(4fi B (B - B )). The black 
curve is for C e = 400^ rec //c rec , where hk rec = ft/A and 
E rec = h 2 k^ ec /2m are the recoil momentum and energy, 
and A = 852 nm is the laser cooling wavelength. In this 
strong-coupling regime, the resonance has a symmetri- 
cal dispersive-like shape. At any given field, all atoms 
within the distribution V r ^ v contribute to it, and the col- 
lision shift reaches the unitarity limit. The green curve 
(C e = 0.4^ rec /^rec) illustrates the weak-coupling regime, 
in which the kinetic energy exceeds the elastic width. In 
this regime, the resonance curve is strongly asymmet- 
ric [24[. For B < Bo, the resonant channel is closed 
and the behavior is similar to the far-detuned strong- 
coupling case. For B > Bo-, the resonant channel is 
open. At a given field, only a fraction of the distribu- 
tion V VyV contributes significantly to the frequency shift 
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because of the narrow elastic width. Consequently, the 
total clock shift is smaller than the unitarity limit value. 
The experimental value C e = Co = 4E rec /k rec (thick red) 
is close to the weak-coupling regime. The resonant be- 
havior of the clock shift is clearly visible, and the inset 
shows that it occurs at the zero-temperature resonant 
field £?o, where this model predicts a singularity even at 
finite temperature. A fit of our model to the measure- 
ments for Mp = 3 (Fig. |4] right) captures the main fea- 
tures of the data, and in particular its asymmetry. This 
fit yields Bo = 5 ± 1 mG. Were the resonance occur- 
ring in the channel, the sign of the clock shift would 
be reversed. Therefore this analysis, independent of our 
coupled-channels results, confirms that the resonance oc- 
curs in the /?7 channel. 

We have measured multiple Feshbach resonances in 
133 Cs at ultralow magnetic fields using a fountain clock, 
and characterized them theoretically using the coupled- 
channels approach. We have identified the resonant 
bound state to be the weakly-bound state of the triplet 
potential and explained their multi-peak structure. They 
have been observed in a regime where the kinetic energy 
dominates over the resonance width, which causes them 
to appear as asymmetric features in the ^-dependence 
of the clock shift, as captured by our finite-temperature 
Monte-Carlo simulations. The resonant triplet state can 
also be brought to resonance using a weak microwave 
field tuned far away from the single-atom resonance [2l| , 
thus leaving the single-atom Physics unaffected, which 
could also be useful for metrological applications. 

We acknowledge many fruitful discussions with J. Dal- 
ibard, P. Rosenbusch, and C. Salomon. 
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